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Abstract 

A major problem in the quantization of fields in curved spacetimes is the ambiguity in the choice of a 
Fock representation for the canonical commutation relations. There exists an infinite number of choices 
leading to different physical predictions. In stationary scenarios, a common strategy is to select a vacuum 
(or a family of unitarily equivalent vacua) by requiring invariance under the spacetime symmetries. When 
stationarity is lost, a natural generalization consists in replacing time invariance by unitarity in the evolution. 
We prove that, when the spatial sections are compact, the criterion of a unitary dynamics, together with the 
invariance under the spatial isometries, suffices to select a unique family of Fock quantizations for a scalar 
field with time dependent mass. 
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In contrast with the situation found in Quantum Mechanics (where one can appeal to the Stone- 
von Neumann uniqueness theorem yj]), linear canonical transformations (even if time indepen- 
dent) are not generally implemented as unitary transformations in Quantum Field Theory (QFT). 
This fact introduces an ambiguity in the choice of a Fock representation for the cannonical com- 
mutation relations (CCRs). Different choices of creation and annihilationlike variables are related 
by linear canonical transformations which cannot be all unitary under quantization. As a conse- 
quence, descriptions which are classically equivalent become inequivalent in the quantum realm. 
The number of linear canonical transformations that cannot be implemented unitarily is infinite, 
actually. Therefore, there exist infinitely many (intrinsically) distinct possibilities for the choice of 
creation and annihilation operators, and hence of inequivalent vacua [2]. This is a severe problem, 
since each of these Fock representations provides a QFT with different physical predictions. 

In order to remove this ambiguity, inherent to QFT in curved spacetimes, the usual strategy is 
to employ the spacetime symmetries of the field theory. Namely, these classical symmetries must 
be inherited in the quantum theory, where they must have a natural unitary implementation. This 
is more stringent than just unitarity: one demands that the structures used in the Fock quantization, 
and hence the vacuum, are invariant under these symmetries. This immediately ensures unitarity, 
while the opposite is not generally true (unitarity does not require vacuum invariance). 

The relevant information on the choice of creation and annihilationlike variables (up to re- 
dundancies like the mixing of only creation variables) is encoded in a basic structure called the 
complex structure. Let us consider the symplectic form of our field system, which encapsules all 
the information about the Poisson brackets. Then, a complex structure is a linear map in phase 
space whose square is minus the identity, leaves invariant the symplectic form and, combined with 
it, provides an inner product in that phase space. Given a complex structure, there is a standard 
procedure (see [2J]) to obtain a Fock representation of the CCR's (or, more precisely, of the corre- 
sponding Weyl relations). Now, we say that a certain group of symmetries admits a natural unitary 
implementation in a Fock quantization if and only if the chosen complex structure is invariant 
under that group. 

The best known example where a unique vacuum is selected by imposing symmetries is QFT in 
Minkowski spacetime: uniqueness follows from the invariance under the Poincare group [2]. But 
in general cases, and in particular in nonstationary settings, there is simply not sufficient symmetry 
to pick up a unique Fock representation. 

In the absence of stationarity, it seems quite natural to minimally relax the requirement of 



invariance under time evolution, replacing it with the condition that the dynamics be unitary. This 
is precisely the criterion, additional to symmetry invariance, that we will put forward in this letter 
to reach uniqueness in the Fock quantization. In fact, this criterion has been successfully tested 



in the context of inhomogeneous cosmologies y 



4] as well as for scalar fields with generic time 



dependent mass on J-spheres, with d = 1,2, 3 I50 (see also |6[] for different criteria concerning 
free fields in 1 + 1 dimensional de Sitter space). 

Specifically, we are going to consider the Fock quantization of a scalar field with generic time 
varying mass (i.e., subject to a time dependent quadratic potential) on a general Riemannian com- 
pact space in three or less (spatial) dimensions. We will prove that a unique, preferred Fock 
representation is selected by imposing the criterion of unitary dynamics and a natural unitary im- 
plementation of the spatial symmetries. 

It is worth commenting that scalar fields propagating in certain nonstationary spacetimes can 
be reformulated, via a time dependent canonical transformation in which the field is scaled by a 
time function, as scalar fields with a time varying mass but in a static background. For instance, 
this occurs in Friedmann-Robertson- Walker spacetimes. Hence, our analysis has immediate ap- 
plications in cosmology, e.g. in the study of cosmological perturbations in models with compact 
spatial topology. This includes the physically important case of flat models with compact sections 
of three-torus topology. 

The description of the scalar field as a Klein-Gordon field with time dependent mass is known 
to have particularly nice properties^ C7Q, and recent results indicate that it can have a privileged 
behavior in terms of its quantum dynamics [8Q. Moreover, it is precisely the field description 
adopted in Mukhanov's formalism for cosmological perturbations [9]. In this context, what we are 
going to show is that, rather than renouncing to a unitary dynamics in cosmology, one can select a 
unique Fock representation precisely by invoking it, together with the spatial symmetries. 

In more detail, we will analyze a real scalar field tp propagating in a globally hyperbolic back- 
groud IxE, where I is a time interval and £ is a Riemannian compact manifold with metric 
h a b. We will study only the case of orthogonal foliations. To pass to the canonical formulation, 
we choose an (arbitrarily) fixed time to, and consider the pairs of data (sp,P) = (<p\ to , Vh(p\ to ), 
where <p denotes the time derivative and h is the determinant of h a b. The phase space is thus the 
1 Unitary dynamics has in fact been achieved in some particular cases, though the uniqueness of the representation 
respecting this unitarity has not been proved before to the best of our knowledge. A detailed discussion of the 
relation between the approach presented here and earlier work on QFT in Friedmann-Robertson- Walker spacetimes 
will be the subject of a future publication. 



set of pairs {(<p,P)}, equipped with the symplectic form Q determined by the Poisson brackets 
{(p(to, x), P(to,y)} = 8{x - y) (independent of the choice of t ). Note that the configuration variable 
(p is a scalar and the momentum P is a scalar density. 

We introduce now the Laplace-Beltrami (LB) operator A associated with the metric h, as well 
as the complex structure: 



in 



-(-hA)- 1 ' 2 
(-M) 1/2 



( \ 



\ j 



(1) 



Clearly, the Fock representation defined by J is the analogue of the free massless field represen- 
tation for Minkowski spacetime. In that case, one can see that J is invariant under the evolution, 
and therefore a natural unitary implementation of the dynamics is achieved then. 

Let us consider the general case of a field with time dependent mass, whose canonical equations 
of motion are 



tp:=-Lp, P= y/h[A<p-f(t)<pl 



(2) 



Here fit) is a rather arbitrary function (apart from some mild conditions specified in [5]). This set 
of equations is equivalent to a Klein-Gordon field equation with square mass equal to f(t). 

We next discuss the necessary and sufficient conditions for a unitary implementation of the 
dynamics dictated by ©, with respect to the Fock representation determined by the complex 
structure Jq. Crucial in this analysis are the properties of the LB operator in compact spaces [10], 
which allow us to perform a mode decomposition of the field, just like with the standard Fourier 
series. In this general setting, the inner product on the space of functions on S is determined by 
the integration with respect to the metric volume element. 

Let then {Y,,/,,} denote a complete set of real orthonormal eigenfunctions of the operator A, 
corresponding to the discrete set of eigenvalues {-oj 2 }, which satisfy oj 2 — * oo as n — » oo. In 
general, there may be some degeneracy, i.e., the eigenspace with eigenvalue (minus) a> 2 may have 
dimension g n greater than the unity (although necessarily finite). This degeneracy is accounted for 
by the label /„ = l,..„g n . In what follows, every sum is performed over the whole spectrum of 
eigenvalues, including degeneracy. 

We then decompose the field (p in terms of eigenmodes: 



<P = J]<lnl n X i'nl,r 



0) 
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The degrees of freedom of the system are thus encoded in a discrete set of real modes q n!n , which 
obey decoupled equations of motion 



q nln + K + f(t)]q„i n = 0. (4) 

Notice that all the modes with the same value of n have the same dynamics, regardless of the 
degeneracy label /„. 

The unitary implementation of the dynamics depends on the behavior of the canonical evolution 
matrices when n — » oo. It is convenient to introduce the annihilationlike variables 

a nln = - — (co n q n i n + ip n ij , (5) 

together with their complex conjugates: the creationlike variables a* n (clearly, this definition is 
meaningless for co n = 0; however, our analysis is not affected by any finite number of modes, 
so we will consider exclusively nonzero modes from now on). Here, p n i n = q n i n is the canonical 
momentum conjugate to q„i n . In these variables, the complex structure J , which has already 
been diagonalized in 2 x 2 blocks by means of the mode decomposition, is further diagonalized 
within each block, taking the standard form Jo(a n i„) = i a ni„ and Jo(ci* n[ ) = ~^ a * n i (this means that 
a„i n and a* nI are precisely the variables quantized as the annihilation and creation operators in the 
corresponding Fock representation). 

The time evolution of the canonical variables (a n i n ,a* tl ), from the fixed reference time t to an 
arbitrary time t, is a linear transformation which is block diagonal, owing to the decoupling of 
modes, and which has the form 

a nln (t) = a n (t, t )a n i n (f ) + J3 n (t, t )a* nln (to), (6) 

given its independence on the degeneracy label. For each label n, the functions a n (t, t ) and fi n (t, to), 
which completely characterize the classical evolution, satisfy 

|ar„(f,/b)l 2 -|ft»(Uo)l 2 = l, VUo- (7) 



As explained in the canonical transformation © is unitarily implementable in the Fock 



representation defined by the complex structure 7o if an d only if 

Yj8n\Pn(t,t )\ 2 <°°. (8) 



Often, this condition is presented as the requirement that the evolved vacuum contain only a finite 
number of particles. 

We can now make use of the asymptotic analysis performed in [5] (see Sec III.B), for the 
harmonic oscillator equation of type ©, in the limit of large co 2 , (i.e., n —* oo). It was shown that, 
for any function /(f), the leading term in fi n (t, f ) is proportional to 1 fo? n . Therefore, condition dU) 
is equivalent to 



Actually, this last condition is satisfied for all Riemannian compact manifolds in three or less 
dimensions. One can see this from well known asymptotic properties of the spectrum of the LB 
operator, which ensure that the number of eigenstates with eigenvalue equal or smaller in norm 



than co 2 - known as the counting function - does not grow faster than a> d [10]. 

Proceeding with our analysis, suppose that the manifold (S, h a b) possesses an isometry group 
G. Then the LB operator is invariant under G, and it follows immediately that this group of trans- 
formations translates into symmetries of the field equations ©, or into canonical transformations 
which commute with the dynamics in the canonical formulation. Consequently, we require those 
symmetries to have a natural unitary implementation in the quantum theory. This is automatically 
granted in the Fock representation defined by the complex structure Jo, since it depends only on 
the metric, and is therefore invariant under the group G. 

Hence, we have a complex structure -Jq- which is invariant under any existing symmetry group 
and provides a representation in which the dynamics dictated by © is unitary. We are now going 
to show that there is no other inequivalent Fock representation with the same properties. Namely, 
we will see that any other complex structure which is invariant under a symmetry group G and 
allows for a unitary implementation of the dynamical evolution, defines a quantum representation 
which is unitarily equivalent to the one fixed by J . 

Let us start with the description of G-invariant complex structures, based on a simple applica- 
tion of Schur's lemma Q3|— 13|] . First, we consider the action of the group G on the set of fields (p, i.e. 
on the configuration space, hereafter called Q. Since G preserves the metric, this action is natu- 
rally unitary (with respect to the inner product on the space of functions <3). Each eigenspace of A 
(corresponding to a given eigenvalue) is itself a representation of G, because G commutes with A. 
Since the action is unitary, a given eigenspace is either irreducible under G or can be decomposed 
into a set of mutually orthogonal irreducible subspaces. Then, in the natural decomposition of Q as 



a direct sum of (finite dimensional) eigenspaces £f" of A, each eigenspace can be decomposed in 
turn as a direct sum of irreducible representations <3" ! " of G. Obviously, the sum of the dimensions 
g„ hi of all such representations Q m " must equal the degeneracy g n , and since the representations are 
at least one-dimensional, we have 1 < g mn < g n . 

We next consider the analogous decompositions of the space V of momentum fields P (here, 
the integration is performed with the inverse volume element). Putting then configuration and 
momentum fields together, we obtain the corresponding decompositions of the phase space Y: 

r = (J) r w " = (J) r m " , r m " = or* e r m " . ( 1 o) 

n n,m„ 

Taking into account that the group G acts in the same way on fields <p and on momentum fields P, 
we easily realize that its action coincides on the irreducible spaces Q m " and V m " . 

As a first application of Schur's lemma II 1211 . we conclude that every G-invariant complex struc- 
ture is block diagonal with respect to the decomposition (flUl) , since invariant transformations do 
not mix different irreducible representations. Thus, an invariant complex structure J must have the 
form 



/ = 0/ m „, (11) 



where each J mn is an invariant complex structure on the corresponding space Y m " . For each of these 
spaces, it is always possible to find a basis formed by configuration variables q nln and correspond- 
ing momenta p nln , associated with a subset of orthonormal eigenmodes as those introduced above 
(recall that the different spaces <3 m " are orthogonal to each other). For a given n, the complete set 
{(Jni„,Pni„} is formed by the union of all such subsets, when all the subspaces Y m " of r w " are taken 
into account. To each J mn , then, there corresponds a matrix characterized by four square blocks: 
Jm„> Jm„i Jm„, an d J^ n - The hivariance conditions on those blocks become invariance conditions 
for matrices in the given irreducible representation of G. Hence, again by Schur's lemma, each of 
the above matrices must be proportional to the identity I, i.e., 

rZ=a m l, J% = b m l, JZ = c m l and J% = d m , I, (12) 

where a mn , b mn , c„ hi , and d mn are real numbers. It follows that invariant complex structures further 
decompose within each sector Y m " into a block-diagonal form such that each block is given by the 
same matrix, namely, a 2-dimensional complex structure which only mixes q n \ n with p n i n for each 
value of /„. 
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To proceed further, it is convenient to switch again from the real basis {q„i n , p„i n } to the complex 
variables a„i n and a* nl , and relate a general invariant complex structure with J Q . One can easily 
show that this relation is always of the form J = KJ K~ l , where K is a symplectic transformation. 
Given our above description of invariant complex structures, the symplectic transformation K can 
also be decomposed in 2 x 2 blocks that are all identical in each space Y m " . So, the information 
necessary to determine a given invariant complex structure is encoded in a discrete set of 2 x 2 
matrices "K mn , one per each irreducible representation of G. They are then totally determined by a 
set {/c m „, AmJ of complex numbers that satisfy the condition |/c m J 2 - \A m J 2 = 1, Vm„ (see e.g. [|5|]). 

On the other hand, it is known that any symplectic transformation R admits a unitary imple- 
mentation with respect to the complex structure J = KJ K~ l if and only if the transformation 
K~ l RK is unitarily implementable with respect to 7o OD. Therefore, the condition for unitary im- 
plementation of the dynamics with respect to a complex structure J = KJ K~ l can be replaced by 
a condition with respect to J . The effect of the transformation K is to replace the functions a n 
and /?„ by new ones, and to possibly lift part of the degeneracy because, in addition to the value 
of n, these new functions can depend on the irreducible representation of G as well. In particular, 
one can compute the relation between the new beta coefficients, hereafter denoted by p J m , and the 
functions a n and J3 n : 

Pt„(t, h) = (K* m ff3 n (t, t ) - A 2 J? n (t, to) + 2iK; n A m 3Mt, t )l (13) 

where the symbol 3 denotes the imaginary part. 

Let us then suppose that, besides J , there exists another invariant complex structure J which 
also allows for a unitary implementation of the dynamics. This unitarity implies that the sequences 
given by y!g^fi J m (t, to) (where we have taken into account the remaining degeneracy) are square 
summable in m n at all possible values of t. We are then in a situation which is completely similar 
to that studied in ||5|], and our result is readily obtained by repeating the arguments presented in 
that work. First, using the fact that condition © is satisfied, it follows from the square summa- 
bility of VgmTAi„(X t ) that the sequences with elements ^[g^,^[a n {t, to)] A m Jk* mn are also square 
summable, V? e I. Moreover, the asymptotic behavior of the sequences {3[a n (t, to)]} was analyzed 
in Q5D (see Sees. IIIB and IVC), and it was shown there that square summability follows as well 
for the sequence 



A m 

— sin 

K m„ 



(O n (t - t ) + 
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(14) 



Finally, we can apply the line of reasoning used in Sec. IVC of [5] -employing Luzin's theo 



rem and integrating over a suitable measurable set in the time interval I- to show that the square 
summability of (fl4l) implies that 2«,m„ gm n \^m n \ 2 < 00 • But it can be easily seen J3, 5] that this 
is precisely the necessary and sufficient condition for the unitary equivalence between the Fock 
representation defined by J and the corresponding one defined by J . This concludes our proof. 

Summarizing, we have shown that the Fock representation determined by J , analog of the 
free massless field representation for Minkowski spacetime, respects the invariance under the spa- 
tial isometries and provides a unitary dynamics for any compact spatial manifold in three or less 
dimensions. Furthermore, among the class of representations defined by invariant complex struc- 
tures, this choice is indeed the unique one (up to unitary equivalence) with a unitary evolution. 

This uniqueness result finds applications in a wide class of systems, including test scalar fields 
with time dependent potentials or, in much more realistic scenarios, (conveniently scaled) cos- 
mological perturbations evolving in a nonstationary, homogeneous spacetime, e.g. in Friedmann- 
Robertson-Walker cosmologies. Our results are valid for perfect fluids with isotropic perturbations 
of the energy-momentum tensor in this class of universes. Furthermore, they can be straightfor- 
wardly extended to tensor perturbations. Besides, preliminary calculations indicate that the pro- 
posed uniqueness criterion is valid even for spinorial fields, which can be treated similarly without 
introducing radically new conceptual difficulties. 

It is worth noting that, in intricate situations where no clear symmetries can be identified in the 
spatial manifold S, one can still resort to unitary groups in the space of square integrable functions 
on 2, constructed from operators that leave invariant the eigenspaces of the LB operator. This more 
general notion of symmetry, however, may force the unitary implementation of transformations 
that arise just from an accidental degeneracy in the spectrum of the LB operator. 

We have addressed the uniqueness problem taking for granted a certain choice of field descrip- 
tion. A complete analysis of the freedom permitted by linear canonical transformations should 
also account for changes of field description obtained by means of time dependent scalings, which 
seem natural to contemplate for fields propagating in nonstationary backgrounds or in cases where 
no fundamental field is fixed from scratch. Concerning this freedom, fortunately, the results of [|80 
strongly support that this ambiguity is also removed by the criterion of symmetry invariance and 
unitary dynamics. 

The proof that the Fock representation determined by J implements the field dynamics in a 
unitary way depends critically on the dimension of the Riemannian manifold. In general, condi- 



tion © is not fulfilled in four or more spatial dimensions (the inverse of the LB operator in the 
orthogonal complement of the zero modes is not Hilbert-Schmidt). In such cases, an open issue 
is whether one can still find a Fock representation which leads to a unitary evolution and analyze 
whether its equivalence class is picked up uniquely by our criterion. Finally, the assumption of 
spatial compactness is essential to avoid infrared problems [2]. However, in cosmology, where 
scales beyond those provided by the Hubble radius are expected to play no significant physical 
role, the infrared behavior can be ignored, and there is no fundamental difference between the 
treatment of compact and noncompact scenarios. In this context, therefore, our uniqueness result 
has complete generality in 3+1 dimensions. 
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